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Abstract 

The sterile neutrino mechanisms for natal neutron star kicks are reanalyzed. It is shown 
that the magnetic field strengths needed for obtaining the observable values of kicks were 
underestimated essentially. Another mechanism with standard neutrinos is discussed where 
the outgoing neutrino flux in a supernova explosion with a strong toroidal magnetic field 
generation causes the field redistribution in "upper" and "lower" hemispheres of the super- 
nova envelope. The resulting magnetic field pressure asymmetry causes the pulsar natal 
kick. 

1 Pulsar proper motion problem 

The problem of large proper velocities of pulsars, born in supernova explosions (pulsar kick), 
has been discussed for more than 40 years. The total list of publications with the observational 
data is very long. We indicate here only the first papers [HE] where the problem was put, 
and the papers where the data were summarized with the samples of 99 pulsars [3] and of 233 
pulsars [3]. Average speed for the sample of 233 pulsars [4j was estimated at the level of 400 
km/s, with more than 15 % having velocities greater than 1000 km/s. The two fastest pulsars 
PSRs B2011+38 and B2224+64 have ~ 1600 km/s. 

It is important that a correlation was finally established between the directions of pulsar 
velocities and of rotation axes. Initially, a conclusion was made in the paper |5j, based on 
an analysis of the set of 29 pulsars, that mechanisms predicting a correlation between the 
rotation axis and the pulsar velocity were ruled out by the observations. However, in the 
paper [gj strong observational evidence was presented for a relationship between the direction 
of a pulsar's motion and its rotation axis. Analysing a set of 25 pulsars which are younger than 
the ones taken in Ref. [5], the authors [6j conclude that 10 pulsars show an offset between the 
velocity vector and the rotation axis, which is either less than lO" or more than 80", a fraction 
that is very unlikely by random chance. 

Obviously, the reason for the initial kick is a kind of an asymmetry in a supernova explosion, 
but the nature of it has not yet disclosed. There were many attempts to explain this asymmetry. 

Numerous attempts to describe the effect in the hydrodynamics of a supernova explosion do 
not explain the large speeds. Three-dimensional simulation of the explosion with the assump- 
tion of initial asymmetry in the supernova core before the collapse, which increases during the 
collapse, leads to the velocity of a pulsar not more than 200 km/s [?]• Multidimensional simu- 
lation by H.-T. Janka e.a. [8] where the explosion anisotropics develop chaotically, resulted in a 
possible pulsar velocity of 10^ km/s. However, there was no correlation between the direction 
of pulsar velocity and the rotation axis direction in this approach. 
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Along with the hydrodynamic approach, there were several different ideas to explain the 
velocities of pulsars, but all of them operated at speeds of scale of 100 km/s: 

i) evolution of close binary systems [U]; 

ii) acceleration of a pulsar within a few months after the explosion due to asymmetric elec- 
tromagnetic radiation caused by the inclination of the magnetic moment with respect to 
the axis of rotation and the displacement of the center of the star [IQj ; 

iii) asymmetric radiation of neutrinos (antineutrinos) in a collapse via the URCA-processes 
in a strong magnetic field of the scale of lO^'* — 10^^ G in a supernova core |llfll3] . 

The neutrino mechanism looks the most interesting. It is known that neutrinos carry away 
99 % of the total supernova energy ~ 3 x 10^^ erg. When the asymmetry is of ~ 3 %, 
neutrinos carry the momentum of ~ 0.03 E/c. The compact explosion remnant, a neutron 
star with a mass ~ IAMq, gets the same momentum. In this case, its velocity can be easily 
estimated as ~ 1000 km/sec. 

However, neutrinos produced in the electroweak processes have small mean free paths in 
matter of the central part of a supernova and may not cause high- velocity pulsars |14fll6j . 

A lively discussion was generated by the idea [17j, under which the neutrino flux asymmetry 
from a protoneutron star arose due to neutrino oscillations in matter and intensive magnetic 
field. The neutrinosphere for Uj- lies inside the neutrinosphere for Ug, and the resonant transition 

— )• jy^ is possible under certain conditions in the region between the neutrinospheres, where 
Ve are entangled in the medium while Ur are "free" to depart. Hence the surface of the resonant 
transition becomes an effective neutrinosphere for Vr- In the presence of a magnetic field, this 
sphere is deformed along the field. Due to the temperature dependence on the radius, the 
anisotropy of the energy flux carried away by neutrinos, arises. This should cause the kick of 
the nascent neutron star. 

The idea of the pulsar kick due to deformed neutrinosphere [T7] raised, however, serious 
criticism [18]: after the neutrinosphere deformation, the surfaces of the constant temperature 
would be deformed also, because just neutrinos provided a thermal equilibrium. And the main 
problem of the model became clear soon: the existence of neutrinos with the mass ~ 100 eV 
was needed. Established restriction on the neutrino mass, mu < 2 eV, "closed" the model. 

There were also attempts to explain large space velocities of young pulsars with using of 
some possible non-standard properties of neutrinos. For example, a mechanism was proposes 
by E. Akhmedov e.a. [19j, of the resonant spin-flavour precession of neutrinos with a transition 
magnetic moment in the magnetic field of a supernova. The asymmetric emission of neutrinos 
was caused by the distortion of the resonance surface due to matter polarisation effects in the 
supernova magnetic field. The requisite values of the field strengths were declared to be of order 
10^® G, and neutrino parameters were taken within the existing experimental bounds. However, 
as was mentioned in the paper [18], in fact the magnetic fields were required in the model [19] 
more than an order of magnitude larger. 

2 The initial pulsar kick and sterile neutrinos 

The sterile neutrinos came on stage in Ref. [20] (see also [21] for details), where the same 
deformed by magnetic field neutrinosphere as in Ref. [17] was discussed, but instead of oscilla- 
tions Vfj^^r ^ i^e, the transitions were considered into "heavy" sterile neutrinos t'^^r o Vs- The 
attractiveness of the model was in an idea that the heavy sterile neutrino (with the mass-scale 
of a few keV) simultaneously solved two problems: providing an initial velocity of pulsars they 
could also play the role of dark matter. 

However, when we have reproduced calculations performed in Refs. |20p21j we have obtained 
that the result for the asymmetry was overvalued in [21j at 15 times. In other words, for 
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the declared asymmetry, the necessary magnetic field strength should be 15 times larger: not 
- 3 X 10^^ G but ~ 4.6 X lO^'^ G. 

Another scenario of using sterile neutrinos for the pulsar kick explanation with off-resonance 
transitions was developed in Ref. p2]. Sterile neutrinos could be born in /3-processes due to 
the neutrino mixing, but the process is suppressed because of smallness of the mixing angle. 
However, they can take a significant amount of energy due to two factors: 

(1) within the core, neutrinos have energies ~ 150 MeV, which is much greater than the 
energy of active neutrinos ~ 20 MeV, emitted from a neutrinosphere; 

(2) the emission occurs from the volume, not from the surface. 

In the presence of a magnetic field, neutrinos are emitted asymmetrically, and this asym- 
metry is maintained because sterile neutrinos are not absorbed, but fly away freely, unlike the 
situation in the approach of Refs. |llfll3j . However, as our analysis shows, the asymmetry was 
overvalued in ^22j at 40 times at least. In other words, for the declared in Ref. [22j asymmetry, 
the magnetic field strength needed should be 40 times larger: not ~ 10^^ G but 4 x 10^'' G. 
In our opinion, the authors [22] made a mistake in calculation of the value ko defined in Eq. (9) 
and presented in Fig. 2 of [22] . By the way, the authors call the value kQ as the fraction of elec- 
trons in the lowest Landau level, while in fact it is the fraction of electron energy squared in the 
lowest Landau level, defining an asymmetry of the neutrino-electron interaction in /3-processes. 
It can be shown that the result of the paper [22] is wrong, both by direct numerical calculation 
and analytically. Really, using Eqs. (9), (10) of Ref. [22] one can transform the expression for 
the value kQ with a good accuracy to the form: 

where B is the magnetic field strength, /^e and T are the chemical potential and the temperature 
of electrons, and Jniv) are the Fermi integrals: 

oo 


Depending on the electron chemical potential and the magnetic field strength, the value ko is 
overestimated in Fig. 2 of |22j by the factor from 40 to 90. 

In the recent e-print [^], G. Kishimoto presented a detailed numerical analysis of the trans- 
formation of active neutrinos to sterile neutrinos through an MSW-like resonance in the pro- 
toneutron star, in order to provide the pulsar kick. However, after correcting a numerical 
mistake in the version 1 of the e-print, it can be seen that the magnetic field strength needed 
for a desirable effect should be taken of the order of 10^^ G. 



3 Back to standard neutrinos? 

A reasonable question arises: if we really need such strong magnetic fields to provide a natal 
neutron star kick with sterile neutrinos, isn't it possible to manage with standard neutrinos? 

As it was already mentioned, an asymmetry of the standard neutrino radiation in a strong 
magnetic field is not a new topic. For example, in the series of papers by our group |24ft26j the 
asymmetry of neutrino emission in a strong magnetic field: 

was analysed, which arised due to parity violation in the neutrino-electron and neutrino-nucleon 
processes. In a strong field of the poloidal type [27H29j . only due to the process u — ?■ i^e~e^ one 
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obtains ^24j: 

where A£ is the characteristic size of the region where the field strength varies insignificantly, 
and E is the neutrino energy averaged over the neutrino spectrum. The asymmetry is seen to 
be not enough to provide the observable neutron star kick with such field strength. 

It should be noted that the mechanism is known of essential enhancement of the mag- 
netic field strength during a supernova explosion. It is the magnetorotational supernova by 
G.S. Bisnovatyi-Kogan |30p31j . a model for generation of the toroidal magnetic field. A poloidal 
magnetic field being enhanced during the supernova core collapse and being frozen in plasma, 
due to the differential rotation, generates a strong toroidal magnetic field which could be in 
order of magnitude greater than the original poloidal field. 

4 Tangential Neutrino Force 

Let us remind what is a possible integral effect of neutrinos on a magnetized plasma. Con- 
sider first the neutrino-electron processes [25]. A complete set of these processes in plasma, 

z^e^ — 7> i^e^, u ve~e^ , ve~e^ ^ , (5) 

lead to the energy and force neutrino fiux impact on plasma: 

{£, ^^) = J\P- P')o. dn. dW , dn. = 1^ ,(^-!!. 1 • 

Here, dW is the total differential probability of all the processes specified in (0), -P and P' 
are the initial and final neutrino four-momenta, the z axis is directed along the magnetic field, 
dnj/ is the initial neutrino density, and T^/ are the effective chemical potential and the 
spectral temperature of the neutrino gas, and the function ^{'d,R) determines the neutrino 
angular distribution, depending on the angle ■!? between the neutrino momentum and the radial 
direction in the star and on the distance R from the center of the star. It should be noted that 
Eq. (I6|) can be used for evaluating the integral effect of neutrinos on plasma in the conditions 
of not very dense plasma, e.g. of a supernova envelope, when an one-interaction approximation 
of a neutrino with plasma is valid. 

Spectral temperatures for different types of neutrinos are estimated to be [32]: 

~ 4 MeV, Tp^ ~ 5 MeV, T^^^^ ~ ~ 8 MeV. (7) 

The probability of the /3 processes (z^e + n -H- e~ + p) is substantially higher than that for 
neutrino-electron processes, so the /3 processes dominate in the energy balance. As a result of 
neutrino heating the plasma, temperature should be very close to the spectral temperature of 
the electron neutrinos, T ~ Tj,^ . 

As it was shown in Refs. [25], the main contributions into the values £ and J^z in & were 
made by /i and r neutrinos and antineutrinos (as a result of the conservation of CP, neutrinos 
and antineutrinos push the plasma in the same direction). This is because in the vicinity of the 
Ve neutrinosphere the spectral temperatures of the other types of neutrinos differ substantially 
from the plasma temperature T ~ T^,^ . 

For numerical estimates we can conveniently express the contribution from i^-e-processes 
with Ue, '^fi,T, i^^i,T into the values £ and J-'z in the following form: 

{£ J J~ z)i>i ~ 



a 



+ c 



(8) 
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where 

(9) 



0.55.10^0^, 



(i) (i) 

The electroweak constants Cy ,C)^ of the effective Lagrangian of the neutrino-electron inter- 
action in Eq. ([8|) are: 



t ^ = 1 + 2 sin2 % , C^) = 1 , C{f = -i + 2 sin^ , 4^'^^ = 



The temperature dependent function has the form: 



oo 







r,-l). (11) 



For electron antineutrinos one obtains ip{1.25) ~ 0.82, while for muon and tau neutrinos and 
antineutrinos the function is ip{2) ~ 38.5. 

A combined effect of all types of neutrinos interacting with electron-positron plasma is: 



Contribution of the neutrino-nucleon processes was evaluated in Refs. [26]. For the param- 
eters of the shell of a supernova: — 0.2, p ~ 10^^^^^ g/cm'^, one obtains (VA^' means both 
urca-processes and i/A^-scattering) 

4^^) c 2.4 X 10^° ( ^ . (13) 

It is important that the contributions of both neutrino-electron and neutrino-nucleon processes 
are of the same sign. The total neutrino force density is: 



J-L*"*"^') ^ 0.6 X 10^1 ( ^ 



Vioiecy 



cm^ 



(14) 



Note that the force density (|14|) is five orders of magnitude lower than the density of the 
gravitational force and thus negligibly influences the radial dynamics of the supernova shell. 
However, when a toroidal magnetic field |30U31j is generated in the shell, the force ()14p directed 
along the field can fairly rapidly (within times of the order of a second □) lead to substantial re- 
distribution of the tangential plasma velocities. Then in two toroids in which the magnetic field 
has opposite directions, the tangential neutrino acceleration of the plasma will have different 
signs relative to the rotational motion of the plasma. This effect can then lead to substantial re- 
distribution of the magnetic field lines, concentrating them predominantly in one of the toroids. 
This leads to considerable asymmetry of the magnetic field energy in the two hemispheres and 
may be responsible for the asymmetric explosion of the supernova which could explain the dis- 
cussed phenomenon of high intrinsic pulsar velocities. In our view it is interesting to model the 
mechanism for toroidal magnetic field generation taking into account the neutrino force action 
on the plasma both via neutrino-nucleon and neutrino-electron processes. 



^We know that the cooling stage of a supernova shell, known as the Kelvin-Helmholz stage, lasts for around 
10 s. 
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5 Neutrino-triggered magnet or ot at ional pulsar natal kick 



Considered neutrino processes in toroidal magnetic field which is frozen in plasma, provide 
the angular acceleration for an element of plasma at a distance R from the rotation axis: 

^ = ^5-i-2xio3i:^(^). (15) 



p R sec 



It means that during the time ~ 1 sec the increase of the angular velocity is 

An - 10^— i-^] . (16) 
sec VlQi^G J ^ ' 

In the one hemisphere, the angular acceleration coincides with the direction of the initial ro- 
tation, while in another hemisphere, they are opposites. A neutrino flux, pushing the plasma, 
torques the toroids in different directions. 

Thus, three stages of a pulsar kick can be identified: 

i) pre-supernova core is collapsing with rotation during 0.1 sec when a strong toroidal mag- 
netic field is generated due to the differential rotation; 

ii) the neutrino outburst, pushing the plasma by the tangential force along the toroidal 
magnetic field which is frozen in plasma, leads to a magnetic field asymmetry: the field 
strength is enhancing in one hemisphere and is decreesing in another one, during ~ 1 sec; 

iii) the pressure difference arising in the two hemispheres, causes the kick to a core. 

According to the momentum conservation, an energetic plasma jet must be formed opposite 
to the pulsar velocity direction H. 

Surely, a detailed multy-dimensional numerical simulation of the process is needed. Let us 
estimate in order of magnitude what to expect. 

A pressure difference arising in the two hemispheres can be evaluated as: 

An ~ — = ^ — , (17) 

where a = 1/137 is the fine-structure constant. The magnetic field pressure causes the plasma 
acceleration: 

. 1.6 X 10-^ f-^)' (HiHe) AO, (18) 

dt sec2Vl0i6Gy V20kmy \ M ) ' ^ ' 



where i?, Q and are the parameters that characterize the region of a strong toroidal magnetic 
field, see Fig. 1. 

Taking for estimation A0 ~ 15^ ~ | , ~ 45" one obtains 

^ ~ 4 X 10^ f-^V f (^-^\ (19) 

dt " sec2Vl0i6Gy' V20kmy' V ^ / 

In fact, the acceleration is not a constant, and an expansion of the magnetic field volume, 
which reduces the magnitude of the field should be taken into account. The magnetic flux 
conservation provides: •pV'^ = const. 

Within the same geometry, one obtains: 

km / Bo \ ( R \ / Az \^/Vl.4M^ ^ 



where Bq is the initial field strength, A z is a distance traveled by a compact remnant of the 
explosion. 



^This remark was made by Hans-Thomas Janka. 



6 




Figure 1: The region of a strong toroidal magnetic field. 



Conclusions 

There are many mechanisms for pulsar natal kick, and the one using sterile neutrinos and 
proposed by A. Kusenko e.a. looks the most attractive. However, as the analysis shows, 
for the declared effect the magnetic field strength should be much larger, not ~ 10^^ G, 
but > 4 X 10^^ G. 

With such strong magnetic fields, it is possible to manage with standard neutrinos. 

As it was shown in the papers by our group, neutrino-electron and neutrino-nucleon 
processes in a strong magnetic field, cause the appearance of a force density acting on 
magnetized plasma along the field direction. 



Tb ^ 0.6 X 10^^ 



f B \ dyne 

Vioi6Gy 



cm^ 



If the strong toroidal magnetic field is generated in the vicinity of the supernova core 
(magnetorotational supernova model by G.S. Bisnovatyi-Kogan) , the neutrino flux, push- 
ing the plasma, torques the toroids in different directions. In the one hemisphere, the 
additional angular acceleration coincides with the direction of the initial rotation, while 
in another hemisphere, they are opposites. 

We stress that just the toroidal magnetic flelds are considered which can be generated in 
order of magnitude greater than the poloidal flelds used in other approaches. 

There arises the magnetic fleld asymmetry in the two hemispheres, and consequently the 
fleld pressure difference providing the pulsar kick acceleration: 

dt ^ sec2Vl0i6Gy V20kmy V M 



and the kick velocity: 



Vkick ^ 600 



km ( Bp \ ( R \ ( Az 



sec VlO^^Gy V20kmy VSkmy V ^ 



/1.4M0^l/2 
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• Because of large acceleration, a pulsar acquires big velocity during very short time, like 
in a shot. One should remember how he operated, being a child, with a cherry-stone after 

eating cherries: pressing it asymmetrically by fingers, he provided a big velocity to that 
compact object. So, we may consider a kind of "Cherry-Stone Shooting" mechanism for 
pulsar natal kick. 

• According to the momentum conservation, an energetic plasma jet must be formed oppo- 
site to the pulsar velocity direction. 

• A detailed multy-dimensional numerical simulation of the process is needed. We would 
believe it should confirm the effect. 

Acknowledgements 

We are grateful to V. A. Rubakov and all the participants of the seminar at the Theoretical 
Physics Department of the Institute for Nuclear Research of the Russian Academy of Sciences 
for useful discussion. We thank H.-T. Janka and A. V. Borisov for interesting remarks and 
C. Kishimoto for correspondence. A. K. expresses his deep gratitude to the organizers of the 
XV-th International School "Particles and Cosmology" for warm hospitality. 

This work was performed in the framework of realization of the Federal Target Program 
"Scientific and Pedagogic Personnel of the Innovation Russia" for 2009 - 2013 (State contract 
no. P2323) and was supported in part by the Ministry of Education and Science of the Russian 
Federation under the Program "Development of the Scientific Potential of the Higher Education" 
(project no. 2.1.1/13011), and by the Russian Foundation for Basic Research (project no. 11- 
02-00394-a). 

References 

[1] I. S. Shklovskii, Astron. Zh. 46, 715 (1969) [Sov. Astron. 13, 562 (1970)] 
[2] J. E. Gunn and J. P. Ostriker, Astrophys. J. 160, 979 (1970). 
[3] A. G. Lyne and D. R. Lorimer, Nature 369, 127 (1994) 

[4] G. Hobbs, D. R. Lorimer, A. G. Lyne and M. Kramer, Mon. Not. R. Astron. Soc. 360, 974 
(2005). 

[5] A. A. Deshpande, R. Ramachandran and V. Radhakrishnan, Astron. Astrophys. 351, 195 
(1999). 

[6] S. Johnston, G. Hobbs, S. Vigeland, M. Kramer, J. M. Weisberg and A. G. Lyne, Mon. 
Not. R. Astron. Soc. 364, 1397 (2005). 

[7] C. L. Fryer, Astrophys. J. 601, L175 (2004). 

[8] L. Seheck, K. Kifonidis, H.-T. Janka and E. Miiller, Astron. Astrophys. 457, 963 (2006). 

[9] J. R. Gott, J. E. Gunn, and J. P. Ostriker, Astrophys. J. Lett. 160, L91 (1970). 

[10] E. R. Harrison and E. Tademaru, Astrophys. J. 201, 447 (1975). 

[11] N. N. Chugai, Pis'ma Astron. Zh. 10, 210 (1984) [Sov. Astron. Lett. 10, 87 (1984)]. 

[12] Yu. M. Loskutov, Pis'ma v ZhETF 39, 438 (1984) [JETP Letters 39, 531 (1984)]; Teor. 
Mat. Fiz. 65, 141 (1985) [Theor. Math. Phys. 65, 1066 (1985)]. 



8 



[13] O. F. Dorofeev, V. N. Rodionov and I. M. Ternov, Pis'ma v ZhETF 40, 159 (1984) [JETP 
Letters 40, 917 (1984)]; Pis'ma Astron. Zh. 11, 302 (1985) [Sov. Astron. Lett. 11, 123 
(1985)]. 

[14] A. Vilenkin, Astrophys. J. 451, 700 (1995). 

[15] D. Lai and Y.-Z. Qian, Astrophys. J. 505, 844 (1998). 

[16] P. Arras and D. Lai, Astrophys. J. 519, 745 (1999). 

[17] A. Kusenko and G. Segre, Phys. Rev. Lett. 77, 4872 (1996). 

[18] H.-T. Janka and G. G. Raffelt, Phys. Rev. D 59, 023005 (1998). 

[19] E. Kh. Akhmedov, A. Lanza and D. W. Sciama, Phys. Rev. D 56, 6117 (1997). 

[20] A. Kusenko and G. Segre, Phys. Lett. B 396, 197 (1997). 

[21] A. Kusenko, Phys. Rep. 481, 1 (2009). 

[22] G. M. Fuller, A. Kusenko, I. Mocioiu and S. Pascoh, Phys. Rev. D 68, 103002 (2003). 
[23] C. Kishimoto. :' arXiv:1101. 13041 

[24] A. V. Kuznetsov and N. V. Mikheev, Phys. Lett. B 394, 123 (1997); Yad. Fiz. 60, 2038 
(1997) [Phys. At. Nucl. 60, 1865 (1997)]. 

[25] A. V. Kuznetsov and N. V. Mikheev, Mod. Phys. Lett. A 14, 2531 (1999); ZhETF 118, 
863 (2000) [JETP 91, 748 (2000)]; in: Particles and Cosmology, Proc. 10th Int. Baksan 
School, eds. E.N. Alexeev e.a., INR, Moscow (2000), p. 44. 

[26] A. A. Gvozdev and I. S. Ognev, Pis'ma v ZhETF 69, 337 (1999) [JETP Letters 69, 365 
(1999)]; ZhETF 121, 1219 (2002) [JETP 94, 1043 (2002)]. 

[27] R. C. Duncan and C. Thompson, Astrophys. J. 392, L9 (1992). 

[28] P. Bocquet, S. Bonazzola, E. Gourgoulhon and J. Novak, Astron. Astrophys. 301, 757 
(1995). 

[29] C. Y. Cardall, M. Prakash and J. M. Lattimer, Astrophys. J. 554, 322 (2001). 

[30] G. S. Bisnovatyi-Kogan, Astron. Zh. 47, 813 (1970) [Sov. Astron. 14, 652 (1971)]. 

[31] N. V. Ardeljan, G. S. Bisnovatyi-Kogan and S. G. Moiseenko, Mon. Not. Roy. Astron. Soc. 
359, 333 (2005). 

[32] V. S. Imshennik and D. K. Nadyozhin, Usp. Fiz. Nauk 156, 561 (1988) [Sov. Sci. Rev., 
Sect. E 8, 1 (1989)]; D. K. Nadyozhin, in: Particles and Cosmology, Proc. 6th Int. Baksan 
School, eds. V.A. Matveev e.a.. World Sci., Singapore (1992), p. 153. 



9 



